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Abstract 

We report a high-statistics measurement of the relative branching fraction B{D^ 
Tr~^Tr~Tr^)/B{D^ — > i^^vr+vr^) using a 532 fb~^ data sample collected with the Belle detector at the 
KEKB asymmetric-energy e~^e~ collider. The measured value of the relative branching fraction is 
^ 7r+7r-7r°)/H(D° ^ K'Tr+Tr^) = (10.12 ± 0.04(stat) ± 0.18(syst)) x 10"^, which has an 
accuracy comparable to the world average. We also present a measurement of the time-integrated 
CP asymmetry in —>■ ir~^ir~Tr^ decay. The result, Acp = (0.43 it 1.30)%, shows no significant 
CP violation. 

PACS numbers: 13.25.Ft, 14.40.Lb 
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INTRODUCTION 

Using a large data sample of decays accumulated with the Belle detector, we obtain 
a precise determination of the —>■ tt+tt^tt^ branching fraction using the — * K^tt^tt'^ 
decay mode for normalization [l|. This study is the first step towards a high-statistics 
Dalitz-plot analysis of the — > ti^h^tt^ decay. The latter could give insight into the 
controversy concerning the S-wave tt+tt^ contribution in D meson decays 0, 13, 0, [sj and 
provide a sensitive measurement of CP violation in the neutral D meson system. Knowl- 
edge of B{D^ -> p'n)/B{D^ K*K), also based on the D'^ -> tt+tt-tt" Dalitz analysis, 
could improve our understanding of the apparent discrepancy of the measured two-body 
branching fractions for D° KK and tttt with theoretical expectations [f?]. A detailed 
study of the — > tt^'k^-k^ decay as well as of other CP-symmetric final states, can 
be used to further improve statistics for the measurement of the angle 03 of the CKM matrix. 

Since both tt^ti~t^^ and K^tx^-k^ involve a neutral pion and the same 

number of charged tracks in the final state, several sources of systematic uncertainties 
nearly cancel in the determination of the relative branching fraction. The method used 
minimizes any dependence on the assumed decay model. The result obtained is compared 
to recent measurements by the CLEO [3| and BaBar collaborations. 

In this study, we also subdivide the same data into and subsamples to calculate 
the time-integrated CP-asymmetry {Aqp) in the Tr+Tr^vr" decay mode. The 

latter study is motivated by the recent measurements of mixing parameters in neutral 
D-meson system j^. The rate of CP violation predicted by the Standard Model reaches 
~ 0.1% in some Cabibbo-suppressed decays such as — >■ 7r"'"7r~7r° (1, |loj. The value 
of Acp in the ti^t[~tt^ decay obtained in the single existing measurement by the 

CLEO collaboration is [5|. We provide a significantly improved measurement 

of Acp{D^ — > 7r+7r^7r°). This measurement is complementary to other measurements of 
Acp in singly-Cabibbo suppressed decay modes (the most sensitive is one by the BaBar 
experiment in K^K~ , Tr"*"-/?" 



EXPERIMENT 

The Belle detector is a large-solid-angle magnetic spectrometer located at the KEKB 
e"'"e~ storage rings, which collide 8.0 GeV electrons with 3.5 GeV positrons to produce 
T(4S) at the energy of 10.58 GeV [12:]. Closest to the interaction point is a silicon vertex 
detector (SVD) surrounded by a 50-layer central drift chamber (CDC), an array of aerogel 
Cherenkov counters (ACC), a barrel- like arrangement of time-of-flight (TOP) scintillation 
counters, and an electromagnetic calorimeter (ECL) comprised of Csl (Tl) crystals. These 
subdetectors are located inside a superconducting solenoid coil that provides a 1.5 T 
magnetic field. An iron flux-return yoke located outside the coil is instrumented to detect 



Kl mesons and identify muons. The detector is described in detail elsewhere |13l . Il4 
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DATA SELECTION 

For this analysis, we used a data sample of 532 fb^^ To reduce backgrounds and also tag 
the flavor of the or decay, we require that the D°'s originate from D* D^n decays. 
A D*^ candidate is reconstructed from a and a low momentum tt where the charge of 
the latter tags the flavor: D*+ ^°7r+g, D*' D^K^g H- D° candidates are 
reconstructed from combinations of two oppositely charged pions (a pion and a kaon in 
the case of — > K~7r^n^) and one neutral pion formed by two photons. In the case of 
multiple candidates, we choose the best candidate using a value based on the vertex 
information of all charged particles, M{D*) — M{D^), and M(7r'') values. A fit in which 
the n^/K^, ir^ momenta are constrained to originate from a common vertex and have the 
nominal mass of the meson is also performed. 

The following kinematic and topological criteria are applied to the charged track candi- 
dates: the distance from the nominal interaction point to the point of closest approach of the 
track is required to be within 0.2 cm in the radial direction and 2.0 cm along the beam direc- 
tion. We also require the transverse momentum of the track to be greater than 0.050 GeV/c, 
to suppress beam background. Kaons and pions are separated by combining the responses 
of the ACC and the TOF with the dE/dx measurement from the CDC to form a likelihood 
C{h), where /i is a pion or a kaon. Charged particles are identified as pions or kaons using 
the likelihood ratio IZ = C{K)/{C{K) + C{7i)). For the identification of a charged pion, 
we require TZ < 0.4; this requirement selects pions with an efficiency of 93% and a kaon 
misidentification probability of 9%. For the identification of charged kaons, the requirement 
is TZ > 0.6; in this case, the efficiency for kaon identification is 86% and the probability to 
misidentify a pion is 4%. We require Ois.hi'iT^ / K^) < 2.2 rad to improve K/ vr separation |l6| , 
where ^lab is the angle between the particle momentum and the z-axis, defined as the direction 
opposite to that of the positron beam. To suppress random combinations of two photons, 
we impose conditions on the energies of the photons constituting the 7r° candidate (£'^(7r°) 
> 0.070 GeV), the two-photon invariant mass (0.120 GeV/c^ < M(77) < 0.150 GeV/c^ 
which corresponds to 2.8 standard deviations (a) in reconstructed M(77)) and the 7r° 
momentum in the laboratory frame {p\ah{'^^) > 0.35GeV/c) to suppress random combi- 
nations of two photons. The mass difference of D* and candidates must satisfy the 
restriction 0.1449 GeV/c^ < M(7rtagVr+7r-7rO) - M(7r+7r-7r°) < 0.1461 GeV/c^ {2a in re- 
constructed M{D*) — M{D^)). The momentum of the D* in the center-of-mass (cms) 
frame of the T(4S) must be in the range 3.0GeV/c < Pcms{D*) < 4.5GeV/c. The lower 
cut is applied to reject D*^s originating from B mesons. The upper cut excludes the re- 
gion of Perns (-D*) with the largest discrepancy between Monte Carlo (MC) simulation and 
data (the difference is taken into account in the systematic error). To eliminate back- 
ground from —>■ K^TT^ — > (7r+7r~)7r° decays, the following veto on M(7r+7r~) is applied: 
0.455 GeV/c^ < M(7r+7r^) < 0.537 GeV/c^ (6.5cr in the reconstructed Ks invariant mass 
resolution). We also require that the tx^-k^-k^ / K^-k^tt^ invariant mass be in the range 1.79- 
1.91 GeV/c^, which corresponds to 5.5cr in the M{D^) resolution. For events passing this 
requirement, the momenta of the final state particles are refitted using the nominal mass 
as a constraint. These refitted momenta are used to calculate Dalitz plot variables as de- 
scribed below. After applying all selection criteria, we find 123.2 x 10^ tt+tt^tt'^ and 
1221.0 X 10^ D° — > fC'vr+vr" events in our data sample. 
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EFFICIENCY CALCULATION 



To obtain reconstruction efficiencies, 22 x 10^ MC events, uniformly distributed over the 
Dalitz plane (DP), were generated for each of the two modes. They were processed using 



the GEANT detector simulation package [17[ and reconstructed with the same selection 



criteria as for data. To take into account the radiative tail in the D'^ invariant mass 
distribution (Fig. 1) due to final state radiation (FSR), the PHOTOS package 3] was used 
for D° — > TT^TC^TT^ /D'^ K'Ti^TT^ at the generator level. 

Differences in the efficiency of particle identification (PID) selection criteria between 
MC and data events are taken into account as correction weights to each signal event. 
They are obtained using a large sample of D* — > D'^vrtag, — > K^t^^ decays, as a function 
of the momentum and polar angle of the decay products. The uncertainties of these 
corrections contribute to the systematic uncertainty of the result. We apply these weights 
only to the kaon in — » K~'k^'k^ and to the corresponding pion (of the same charge) in 

— > 7r"'"7r~7r°, since the corrections to the PID efficiency for the remaining decay pion and 
the tagging pion (vrtag) cancel in the ratio. 

A certain portion of signal MC events (~ 15%) is reconstructed with one or more 
random particles (7's, 7r°'s, tt^ or K^) combined with true signal daughters. We dis- 
tinguish correctly reconstructed and misreconstructed signal MC events by comparing 
the reconstructed momenta of all the final state particles to the corresponding generator 
information. The correctly reconstructed MC events are used to calculate the reconstruction 
efficiency, and the misreconstructed decays are treated as an additional source of background. 

The M{D^) distributions for correctly reconstructed signal MC events (Fig. la) are fitted 
with a double hyperbolic Gaussian jioj and one regular Gaussian. The M(D°) distributions 
for the misreconstructed signal MC events (Fig. lb) are fitted with a triple Gaussian. The 
results of the M{D^) ffis are used to fix the shape of the signal and misreconstructed signal 
events for the data M(D°) ffi. 

The fraction of correctly reconstructed events in a certain bin depends on its position on 
the DP, i.e. M^(/i7r) vs. M^(7r7r°) {h is or vr). To determine the reconstruction efficiency 
we divide the DP into bins of size 0.1 GeV^/c^ x 0.1 GeV^/c^, accumulate signal MC events 
from the M{D^) signal region, and then normalize by the number of generated events in 
each bin. The calculated values are later used as reciprocal weights for the corresponding 
data distribution. This method takes into account variations of the DP data density and 
minimizes decay model dependence. 



BACKGROUND STUDY 

To describe the shape of background in the M{D^) signal region for —>■ tt+tt^tt^ 
and K'tt'^tt'^, a sample of generic MC events (including all significant processes 

in e^e^ production of T(4S'), uu, dd, ss and cc at the given -y/i), equivalent to ~600 
fb~^, was processed with the same selection criteria as data. All generic MC events recon- 
structed as 7r+7r^7r° were separated into three types: contributions from uu, dd, and 
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FIG. 1: (a) correctly reconstructed (x^/n.d./. = 1.7) and (b) misreconstructed (x^/n.d./. = 1.0) 
signal MC distributions. 

ss fragmentation, and T(4S') BB events; a contribution from D* D°(fC~7r"*"7r°)7rtag 
(cc events) where the charged kaon is misidentified as a pion (the largest source of back- 
ground); and a contribution from cc background that does not involve particle misiden- 
tification and from which the signal is excluded (Fig. 2a - 2c). For the K~7r'^7r^ 
case, there are contributions from uu, dd, ss and T{4S) BB events, a contribution 
from D* — > D''(fC^7r+7r°7r°)7rtag via K*p and — > aiK, as well as a small residual 

e"'"e~ — > cc background (Fig. 2d - 2f). The small peak in the signal region of the latter 
(Fig. 2f) is mainly due to combinations of a D*^ and a random TTtag and has to be taken 
into account. This background is also present in — > 7r'^7r~7r°. As described previously, 
the contributions of misreconstructed signal MC events are treated as separate sources of 
background for both decay modes. 



DATA M(D0) fit 

The M(D°) distribution in data is fitted using fixed MC shapes for the various background 
components, and a signal shape that allows for data-MC differences. The fit function for 
7r+7r-7r° is 

-^1 — -^sig ^ -Psig('-'^add) ~\~ -/Vmisrec ^ -fmisrec ~l~ ^^\^ 

+ Nudsb X Pudsb + -^misid X -Pmisid + Nc X Pc, 

where Pgig and Pmisrec are the shapes of the M{D^) distributions for correctly reconstructed 
and misreconstructed signal MC events obtained from the corresponding MC distributions. 

The M(D°) shapes of m, s-quark and BB decays, misidentified K~7r~^7i^ decays and 

other c-quark contributions are denoted as Pudsb, -Pmisid and Pc, respectively; A^sigj ^misrec, 
^udsb, -^misid and Nc are the normalizations of all the event types and are free parameters in 
the fit. The additional free parameter Ax represents a common shift in the central value of 
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FIG. 2: M{D^) distributions for MC background events in (a-c) the 7r+7r~7r° and (d-f) 

the K~TT~^7T^ sample: (a) K~Tr~^Tr^ events, (b) other e+e^ cc contributions, 

(c) contributions from hght quark and BB decays, (d) a^K K^tt^'k^it^ events, (e) 

K*p K~7r~^Tr^TT^ events, (f) other e~^e~ cc contributions and contributions from light 
quark and BB decays. Events from the M{D^) signal region (1.79 to 1.91 GeV) are selected for 
the branching fraction calculation. 

the Gaussians describing the signal. Similarly, a^dd is a free parameter added in quadrature 
to all the widths of the Gaussian functions (0.3 MeV for M(L)° 7r+7r-7r°) and 0.1 MeV 
for M(D° K-T^+TT^)). The M{D^ K-t:+tt^) fit function has a similar form: 

-^2 — -^sig ^ -Psig('-'^add) ~l~ -^misrec ^ -fmisrec ~l~ /„\ 

Pudsb 

where Pk*p and Pa^K are the shapes of the contributions of the — > K*p and — > aiK 
decays to the M{D^) distribution, respectively, and Nk*p, N^-^k are their floating nor- 
malizations. All other variables are the same as in Eq. (1). Figure 3a shows the fit for 
M{D^ iT~^'iT~7r^) described above, while Figure 3b shows that for M(L>° K' 
The low fit quality {x^/n.d.f. = 3.2) of the latter is due to the large statistics of the 
signal as well as the poor agreement of data and MG simulation for the D° — > Kmr^TT^ 
background. This discrepancy is taken into account as a systematic error due to the fit 
uncertainty. 



CALCULATION OF THE SIGNAL YIELDS 

After the parameters of the data M{D^) distributions are obtained from the fit, we fill 
separate M^(/i7r) vs. M^(7r7r°) Dalitz histograms with events from the signal M{D^) region 
for data and simulated background with the normalizations fixed from the fit (Fig. 4). The 
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FIG. 3: (a) M{D^ 7r+7r~7r°) data fit {x^/n.d.f. = 1.5). Data is represented by the points and 

the curve is the fitted sum of all the contribTitions (simulated signal and background) . The vertical 
dashed lines indicate the M{D^) signal region. Background: misreconstructcd signal (dashed line), 
"7r"'"7r'' with a misidentified kaon (shaded histogram) and other sources, i.e. other cc and 
light quark contributions (hatched histogram), (b) M{D^ — K'-k^-k^) data fit [x^ /n.d.f. = 3.2). 
Background: misreconstructed signal (dashed) and other sources (hatched). The fit results shown 
correspond to the second step of the B{D^ -k^ ir^ ir^) / B{D^ — K'tt'^it^) calculation, which 
takes into account the decay model (see text). 



bin size is 0.1 GeV^/c'^ x 0.1 GeV^/c^, the same as for signal MC events. 

The number of — > 7r'*"7r~7r° signal events in each bin is calculated as follows: 

— ly -iVmisrec ^ '-'misrcc ^g-^ 

- ^udsh X B\dsh - -^misid X -Smisid " -^c X -B*, 

where is the number of data events, and A^misrec x S\^^^^^ ^udsh x B\idsh-> x Bl 
are the numbers of different background events in the i-th bin. The procedure for the 
D° K'TT^TT^ case is similar. 



The total number of signal events for both decays is obtained by summing the number 
y of events over all bins: S = ^y/e*, where the reconstruction efficiency in each bin 

= ^^icc/^gcn used as a reciprocal weight {nl^^ and rigg^ are the numbers of reconstructed 
and generated events in the i-th bin). 

At this point, we return to the stage of obtaining the M(D^) distributions from signal 

MC simulation and perform another iteration of the same procedure using the Dalitz 
histogram for data as an approximation of the D° decay model (for each of the two decay 
modes). As mentioned above, signal MC events, used at the first step of our calculations 
(the entire procedure described above), are distributed uniformly over the DP. At the 
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FIG. 4: D° 7r+7r-7r°: Dalitz {M^{Tr+TT-) vs. M2(7r+7r°)) distributions for (a) efficiency, (b) 
simulated background and (c) data. 

second step, the decay model is taken into account by weighting entries in a histogram 
according to their positions on the DP. This is done to obtain a more exact M{D^) 
distribution for the signal and misreconstructed signal MC events. The distributions are 
refitted and the resulting background normalizations in Eq. (2) are used to recalculate the 
signal yields Y' in Eq. (3). This results in S{D^ 7r+7r-7r°) = (2403.6 ± 9.2) x 10^ and 
S{D^ K-TT+TT^) = (23751 ± 24) x 10^. 



SYSTEMATIC UNCERTAINTIES 

The sources of systematic uncertainty are as follows; the values quoted are relative 
fractions. The estimate of the error due to the tracking efficiency uncertainty is based 
on a large sample of partially reconstructed D* — > D'^TTtag, -D° fC^Tr^Tr" decays. The 
uncertainty for the two charged tracks — tt+tt" or K~ti~^ — cancels to a large extent in 
the ratio of the D° — > 7r"'"7r~7r° and — > K'n'^n^ branching fractions. It contributes only 
0.01% to the overall systematic uncertainty. We assume that the 7r° and the tagging pion 
(from the D*) reconstruction efficiencies fully cancel in the ratio of the branching fractions. 

The uncertainties of the corrections to the efficiency of PID selection criteria contribute 
±0.91% to the systematic uncertainty of the result. The statistical error of the signal 
MC sample contributes ±0.30% to the total systematic uncertainty The systematic 
uncertainty due to the fractions of signal and various backgrounds, which are fixed from 
the M{D^ — i> 7r+7r^7r°) fit results, was determined by varying the fractions within their 
errors (±0.61%). The correlations between the fit parameters were accounted for using the 
covariance matrix obtained from the fit. The uncertainty due to the M(D° — > K~'k^t:^) 
fit (±0.30%) was estimated by relaxing or fixing relative normalizations of some of the 
background types. 

Our method for calculating B{D^ — > 7r'^7r''7r^) /B{D^ — > K'n'^n^) minimizes the 
uncertainty due to modelling the n~^n~7r'^ and K~ti^ti^ decays. However, 

the model dependence of the background is included in the total systematics. The 
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level of background in the K~7r~^7T^ decay is small and its effect on the ratio of 

— > TT'^TT~7r^)/B{D'^ — > K^-K^TT^) IS negligible. The dominant source of background 
for the — * tt+tt^tt^ mode is the — ^ K'n'^'n'^ decay with the kaon misidentified as 
a pion. The normalizations of the K^ti^ti^ submodes {D^ pK, K*it, 

K*°n° and non-resonant K tt+tt'^) are varied within the uncertainties 

of their branching fractions [20| and the resulting differences from the central value 
of B{D^ — > 7r+7r~7r°)/i3(D° — >• K~n^TT^), summed in quadrature, are treated as the 
background model uncertainty (±0.48%). 

Changing the DP bin size from 0.1 GeVVc^ x 0.1 GeVVc^ to 0.05 GeVVc^ x 0.05 
GeVVc^ yields a ±0.54% difference in the value of B{D^ n+n~7r^)/B{D^ K-jt+tc'^). 
We study the effect of the selection criteria upon the fraction of correctly reconstructed 
signal MC events and obtain a corresponding error of ±0.10%. We varied the event selection 
criteria in order to estimate the systematic error due to any inadequacies in the background 
description. Varying the Ks veto yields a ±0.50% systematic uncertainty. The change of 
B{D^ 7r+7r~7r°)/i3(L)° K~7T~^7r^) due to the variation of the Pcms{D*) upper cut is 
neghgible. Varying the Pcms(-D*) lower cut yields a relatively large uncertainty of ±0.77%. 
Uncertainties due to the variation of other selection requirements are listed in Table [H The 
total uncertainty is obtained by adding all contributions in quadrature. 

TABLE I: Contributions to the relative systematic error on B{D'^ tt~^tt~it^)/B{D^ i^~7r+7r°) 



Source 


Error, % 


Source 


Error, % 


PID corrections 


0.91 


Selection criteria: 




MC statistics 


0.30 


Ks veto 


0.50 


Fitpo ^ 7r+7r-7r°) 


0.61 


PcmsiD*) 


0.77 




0.30 


M(K-7r+7r737r) 


0.36 


TT^TT^ir^ backgr. model 


0.48 


AM 


0.30 


Binning 


0.54 




0.40 


MC misreconstruction 


0.10 


M(7r0) 


0.20 


Tracking 


0.01 


Plab(7r°) 


0.16 


Total 


1.79 



RESULTS FOR BiD° tt+tt- tt^) /B{D^ R-tt+tt^) 

Summarizing the discussion above, we obtain the following ratio of the branching frac- 
tions: 

B{D^ K-TT+n^) ~ S{DO R-tt+tt^) ~ (4) 
= 0.10120 ± 0.00040(stat) ±0.00181(syst) = 0.1012 ± 0.0019. 

We can compare our measurement of the ratio with a recent result obtained by BaBar . 
There is a 2a difference between the central values; the accuracies of the measurements 
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are comparable. To compare results from different experiments, we multiply the obtained 
value of Eq. (4) by the 2007 world average of ^ R-n+Ti'^) = (13.5 ± 0.6)% 



to calculate the absolute branching fraction for the — >■ 7r+7r^7r° decay (Table [IT])- In 
a recent study by CLEO 0], the relative branching fraction B{D^ -k^ -k^ -k^) / B{D^ — >■ 
K^7r+) is measured to be 0.344 ± 0.005(stat) ± 0.012(syst). Using the world average value 
of B{D^ — > K~Ti^) = (3.82 ± 0.07)% from j20|, one can calculate the absolute branching 
fraction of B{D^ 7r+7r~7r°) from CLEO data as shown in Table [TTl A comparison of the 
corresponding values for the absolute branching fraction B{D^ —>■ 7r+7r~7r°) shows that the 
results are in good agreement (2l| . 



TABLE II: ^3(1?° vr+Tr-vr") by Belle, BaBar @] and CLEO 0]. The first two errors are statistical 
and systematic, respectively, and the third one (the fourth column) is the normalization uncertainty. 
The latter is common in the Belle and BaBar results. 



Group 


Nev, 10^ 


B{D° ^ ^+7r-7r0)/e(D0 ^ R-tt+tt^) 


^(L»o^7r+7r-7rO),10-3 


Belle 

BaBar 

CLEO 


123.19 ±0.49 
60.43 ± 0.34 
10.83 ±0.16 


0.1012 ±0.0004 ±0.0018 
0.1059 ±0.0006 ±0.0013 


13.66 ±0.05 ±0.24 ±0.61 
14.30 ±0.08 ±0.18 ±0.64 
13.14 ±0.19 ±0.46 ±0.24 



MEASUREMENT OF Acp 

We subdivide the n'^n^n^ sample into n^n^n^ and D° tt+tt^tt^ subsamples 

to calculate the value of the time- integrated CP-asymmetry using the same method for 
calculating the signal yield that was used for the relative branching fraction. The fitted 
M{D^) distributions of the data are shown in Fig. 5. The resulting values of S* = ^yYe* 
are 

Sdo = (1154.7 ± 6.7) x 10^ 

S^o = (1144.7 ± 6.6) X lOl ^ ^ 

Their sum differs from the value used to calculate the branching fraction, because the 
corrections for the PID efficiency of the pions originating from the cancel out in the 
case of the Acp calculation and thus are not applied. 

A detector bias may exist that leads to different efficiencies for reconstructing positively 
and negatively charged tracks. This may be due to charge-dependent effects such as 
opposite signs of the Lorentz angle with respect to the curvature of tracks in the CDC, 
and a difference in nuclear interactions with the detector material for positive and negative 
tracks. The former is partially taken into account by generating signal MC samples for 
— i> D^TT^^g and D*~ D^vr^g separately. However, nuclear interactions of charged 
tracks with the detector material are imperfectly simulated. This fact also causes a 
systematic difference between tracking efficiencies for positive and negative particles and 
has to be taken into account. 

Since D'^ and are distinguished only by vr^^g, and the neutral D meson decays are 
charge and particle-type balanced, the uncertainties of the reconstruction efficiencies of the 
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FIG. 5: (a) M{D^ Tr+vr^vr'^) and (b) M{D^ vr+Tr^vr'^) data. Background: misreconstructed 
signal (dashed line), — > K^ti^tt^ with misidentified kaon (shaded histogram) and other sources 
(hatched histogram). Events from the M{D^) signal region (1.79 to 1.91 GeV/c^) are selected. 

pions originating from do not affect the result. We consider the uncertainties in the track- 
ing and PID efficiencies of the tagging pions as the main source of systematic errors for Acp- 

The uncertainty of the tracking efficiency was obtained using the same method used for 
the systematics of the —>■ TT~^7T~7r^ /D^ K~7i~^7t^ ratio, but in this case, positive and 
negative vrtag's were treated separately. The calculation of the systematic error takes into 
account the momentum dependence. The errors for vr^t^g and n^^^ were propagated to Acp 
assuming them to be uncorrelated. The charge- dependent data/MC PID corrections for 
TTtag were obtained using independent D* D^{Ks'n'^)7!'tag data and MC samples. 

In general, the D-meson distribution is an asymmetric function of cos(^) (where 6 is 
the polar angle) due to the interference of virtual 7 and in the process of c-quark 
pair production. If the detector acceptance in the center-of-mass frame were perfectly 
symmetric, the cos(6') dependent asymmetry of and (D~^ and D~ etc.) production 
would cancel out in the integral over cos(6') in a symmetric interval. However, the detector 
acceptance is not symmetric and a possible forward-backward asymmetry {Aft,) should 
be taken into account. A data sample of K^K^ ,tt^7i~ decay events was used 

to calculate Afb[cos{6)). This function was then used to weight the MC tt+tt^tt^ 
distribution, which was then normalized to the total number of MC —>■ vr'^Tr^Tr" events. 
The calculated value equals 0.15% and is treated as the systematic uncertainty related to 
the forward-backward asymmetry. Other individual sources of systematic uncertainties are 
listed in Table IIIII Systematic errors for each flavor are calculated similarly to those for 
i5(D° — ^ 7r"'"7r~7r°)/i3(D° K~7i~^7i^), propagated to Acp, and then added in quadrature. 
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TABLE III: Systematic uncertainties for Acp: 



Source 


MC Stat. 


Tracking 


Fit 


Ks veto 


PID 


Binning 




Total 


a, % 


0.24 


1.01 


0.58 


0.23 


0.15 


0.05 


0.15 


1.23 



The resulting value of the asymmetry is 

^CP = (Soo - S^o)/{Sdo + S^o) = , , 

= (0.43±0.41(stat) ± l.Ol(track) ±0.70(other syst))% = (0.43 ± 1.30)%. ^ ^ 

This result is consistent with CP conservation in this decay mode; its sensitivity is a signif- 
icant improvement over that of the previous measurement, {lt.g^)% 



SUMMARY 

Using 532 fb~^ of data collected with the Belle detector, a high-precision mea- 
surement of the relative branching fraction B{D^ —>■ TT^7T~7r^)/B{D^ K~7t~^tt^) = 
0.1012 ± 0.0004 ± 0.0018 has been performed. The method applied minimizes possible 
systematic uncertainties due to the decay model. The mode —>■ K'tt^h^ is 
chosen for normalization to avoid most of the tracking and particle identification un- 
certainties. We also calculate the value of the time-integrated CP asymmetry to be 
Acp{D^ — > 7r+7r~7r'') = (0.43 ± 1.30)%, which is consistent with zero. The sensitivity is 
significantly better than that of the previous measurement 



ACKNOWLEDGEMENTS 

We thank the KEKB group for the excellent operation of the accelerator, the KEK 
cryogenics group for the efficient operation of the solenoid, and the KEK computer group 
and the National Institute of Informatics for valuable computing and Super-SINET network 
support. We acknowledge support from the Ministry of Education, Culture, Sports, Science, 
and Technology of Japan and the Japan Society for the Promotion of Science; the Australian 
Research Council and the Australian Department of Education, Science and Training; the 
National Natural Science Foundation of China under contract No. 10575109 and 10775142; 
the Department of Science and Technology of India; the BK21 program of the Ministry of 
Education of Korea, the CHEP SRC program and Basic Research program (grant No. ROl- 
2005-000-10089-0) of the Korea Science and Engineering Foundation, and the Pure Basic 
Research Group program of the Korea Research Foundation; the Polish State Committee 
for Scientific Research; the Ministry of Education and Science of the Russian Federation and 
the Russian Federal Agency for Atomic Energy; the Slovenian Research Agency; the Swiss 
National Science Foundation; the National Science Council and the Ministry of Education 
of Taiwan; and the U.S. Department of Energy. 



14 



[1] Unless specified otherwise, both flavors of mesons are implied: — > K tt~^it^ and 

K+TT-TT^. 

[2] E.M. Aitala et al. (E791 Collaboration), Phys. Rev. Lett. 86, 770 (2001). 
[3] J.M. Link et al. (FOCUS Collaboration), Phys. Lett. B 585, 200 (2004). 
[4] H. Muramatsu et al, Phys. Rev. Lett. 89, 251802 (2002). 

[5] D. Cronin-Hennessy et al. (CLEO Collaboration), Phys. Rev. D 72, 031102 (2005). 

[6] I.I. Bigi, arXiv:hep-ph/0703132v2, and references therein. 

[7] P. Rubin et al. (CLEO Collaboration), Phys. Rev. Lett. 96, 081802 (2006). 

[8] B. Aubert et al. (BaBar Cohaboration) , Phys. Rev. D 74, 091102 (2006). 

[9] M. Staric et al. (Belle Collaboration), Phys. Rev. Lett. 98, 211803 (2007); L.M. Zhang et al. 
(Belle Cohaboration), Phys. Rev. Lett. 99, 131803 (2007); B. Aubert et al. (BaBar Collabora- 
tion), Phys. Rev. Lett. 98, 211802 (2007); B. Aubert et al. (BaBaR Collaboration), arXiv:hep- 
ex/0712.2249vl; T. Aaltonen et al. (CDF Collaboration), arXiv:hep-ex/0712.1567vl. 
[10] Y. Grossman, A.L. Kagan, Y. Mr, Phys. Rev. D 75, 036008 (2007). 
[11] B. Aubert et al. (BaBar Collaboration), arXiv:hep-ex/0709.2715vl. 

[12] S. Kurokawa and E. Kikutani, Nucl. Instr. and Meth. A 499, 1 (2003), and references therein. 

[13] A. Abashian et al. (Belle Collaboration), Nucl. Instr. and Meth. A 479, 117 (2002). 

[14] Z. Natkaniec et al. (Belle SVD2 Group), Nucl. Instr. and Meth. A 560, 1 (2006). 

[15] Z)*'s originate mainly from continuum. Although we do not apply any topological cuts, the 

yield of D*'s coming from B mesons is negligible: they are rejected by kinematic cuts, mainly 

by the stringent pcms{D*) requirement. 
[16] The backward end of the ACC corresponds to ^lab ~ 2.2 rad, and as a result, kaon-pion 

separation is less efficient for ^lab > 2.2 rad. 
[17] GEANT, R. Brun et al., GEANT 3.21, CERN Report DD/EE/84-1, 1984. 
[18] P. Golonka and Z. W§s, Eur. Phys. J. C 50, 53 (2007). 

[19] / = e^^y^°'^ ■ "K^ ^o)^ gJ_^ where y = \J o? ^ a - (Jyx — XQ)/a)'^. The function tends to a 
Gaussian for |x — xq] ^ \foL ■ a and has exponential behaviour for — xo| ^ \/a • ex. 

[20] W.-M. Yao et al. (Particle Data Group), J. Phys. G 33, 1 (2006) and 2007 partial update for 
edition 2008. 

[21] Our conclusions do not change if instead of the world average value of B{D^ K~7r^) we use 
the value of this branching fraction from the recent high precision measurement of CLEO 2M ■ 
[22] S. Dobbs et al, Phys. Rev. D 76, 112001 (2007). 



15 



